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Abstract 

We describe a method of control of chaos that occurs in area-preserving maps. This 
method is based on small modifications of the original map by addition of a small 
control term. We apply this control technique to the standard map and to the 
tokamap. 
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1 Introduction 



Chaos that arises naturally in Hamiltonian systems with mixed phase space 
is a source of unstable behaviors when resonances overlap. Achieving the con- 
trol of these systems by reducing or even suppressing chaos is a long standing 
and crucial problem in many branches of physics. Among the various methods 
which have been designed to control Hamiltonian chaos, two different frame- 
works must be distinguished : Hamiltonian flows (continuous time) [1,2], and 
symplectic maps (discrete time) [3,4,5,6,7,8,9,10,11,12]. 

A method for controlling Hamiltonian flows has been proposed in Refs. [13,14,15,16]. 
Its aim was to find a small control term / for the perturbed Hamiltonian 
H = Hq + eV where H is integrable, in order to have a more regular dy- 
namics for the controlled Hamiltonian H + eV + /. An explicit algorithm 
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to construct a control term of order e 2 making the controlled Hamiltonian 
integrable was devised. The main advantage of this approach to control is the 
robustness of the control term /, which means that one can tailor another 
control term / which is close to / such that (i) the controlled Hamiltonian 
H + eV + / is still more regular than H + eV and (ii) this control terms 
satisfies some specific requirements of the experiment like a partial knowledge 
of the potential V [14] or the localization of the control in a phase space re- 
gion [16]. Let us stress that this method differs from other methods by the 
fact that the controlled dynamics is Hamiltonian : This makes it relevant to 
the control of inherently Hamiltonian systems [17] such as beams of charged 
test particles in electrostatic waves [18,19], two-dimensional Euler flows [20] 
or the geometry of magnetic field lines [21,22]. 

In this article, the problem we address is to set up a control theory in a sim- 
ilar way for symplectic maps. We consider two-dimensional symplectic maps 
(A, tp) i— > (A', tp') = F(A, tp) on the cylinder KxT which are e-close to inte- 
grability. Our aim is to find a small control term / such that the controlled 
symplectic map F + f is more regular (closer to integrability) than F. In what 
follows, we develop algorithms for finding control terms of order e 2 . 

We describe two approaches to the control of symplectic maps : One approach 
(in section 2) is designed for maps which are generated by Lie transforms. We 
show that it is possible to obtain an explicit expression for a formal control 
term which makes the controlled map integrable. This approach extends in 
an obvious way to 2iV-dimensional symplectic maps. However, since this ap- 
proach is difficult to handle in practice, we develop a second one (in section 3) 
for maps which are written in mixed coordinates (i.e. half implicitly, half ex- 
plicitly). This second method is used for practical purposes since numerous 
area-preserving maps are written in such form. 

In section 4 we apply our control method to two examples : the standard map 
to benchmark our method, and the tokamap, which imposes an additional 
restriction on the type of control. We compute analytically the control term 
and show numerically the phase portrait of these maps and the associated 
controlled maps. As the formal expansion may generate undesirable side-effects 
in phase space, we also discuss suitable modifications of the control terms. 



2 Maps expressed as Lie transforms 

Let A be the vector space of C°° real functions defined on the phase space 
KxT. For H E A, let {H} be the linear operator acting on A such that 

{H}H' = {H,H>}, 



2 



for any H' G A, where {• , •} is the Poisson bracket. Hence A is a Lie alge- 
bra. For two-dimensional maps written in action-angle coordinates, the linear 
operator {H} writes 

OH d OH d 

We consider area-preserving maps generated by a Lie transform 

e fH »e fyi | A | , (1) 

v . 

where H and V are two elements of the algebra A and the linear operator 
e^ is defined by the expansion 
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{H} = V- { H Y 

n=0 U - 



In what follows, we will use two warped additions, one for linear operators 
acting on A (elements of L(A)) and another one for elements of A : First, we 
define the warped addition © between two elements a and j3 of L(A) as 

An explicit series for a © (3 is given by the Baker-Campbell-Hausdorff for- 
mula [23]. Specializing to a = {H} and f3 = {V} for H,V G A, one can 
see that {H} © {V} is also a Lie operator : For instance, the first terms of 
{H} © {V} are 



{H} © {V} = {H} + {V} + \ ({H}{V} - {V}{H}) + -.., 

= {H + V+ l -{H}V + ■■■}, 

where the second equation follows from the Jacobi identity. We associate with 
{H} © {V} an element of A which we denote H EH V if it exists (it does, e.g., 
when H and V are analytic in A and ip, and V small enough) : 

{H}®{V} = {HmV}, (3) 

or equivalently 

e {Hmv} = e {F} e {n ) 

where the first terms of the expansion of H EB V are 

Hmv = H + V + -{H}V + ■■■ . 

Both warped additions © and EB are associative but not commutative. Also, 
we have the following two properties : W = H EB V implies H = W EB (— V), 
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and (-H) ffl (-V) = -{V ffl H). The map given by Eq. (1) will be called the 
map generated by H EE V. 

For instance, if H depends only on the action A, and V only on the angle ip, 
the map obtained from Eq. (1) is the composition of the following two maps 




-d v V{<p)d A 




f A-d v V{<p)\ 



J 



and 




A' 

tp' + u(A') 



where u(A) = 8aH and &a and d v denote the partial derivatives with respect 
to A and (p. Thus 



A" = A-d v V{<p), 

ip" = ip + u;(A-d^V(ip)). 

When V depends also on the action A, the map is much more difficult to write 
explicitly. 

From the operator {H}, we define the linear operator 



H = 1 



This operator 7i is not invertible : for instance TiH = 0. Hence we consider a 
pseudo-inverse of H (if it exists), i.e. a linear operator r on A such that 



H 2 T = n. 

From this operator, we define the operators 1Z and J\f as 

M = HT, 
TZ—l—Af. 



(4) 



Note that H1Z — 0. If T commutes with H, the operators 1Z and Af are 
projectors. 

From now on, we assume that H depends only on the action A. Then the 
pseudo-inverse Y exists formally. The operators {H}, H, T, 1Z and Af act on 
a function U(A, tp) = £ fcGZ U k {A)e ikLp as 
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{H}U = J2*ktu(A)U k (A)e^, 

k 

HU = 2(1 - e- MA)k )U k {A)e tk *, 
k 

ru- V u k{A) ikv 

u i_ e -iuj(A)k K ' 

k S.t. w(A)fc^27rZ 

7H7 = ^ ^(A)e^, 

fe S.t. w(A)feG27rZ 

A/77 = ^ ^(^)e^. 

fc S.t. w(A)fc£27rZ 

Proposition 1 There exists a control term f of order V 2 such that the map 
generated by H ffl V ffl / is canonically conjugate to the map generated by 
H ffl 7£V The control term is given by the expression 

f = (—V) ffl (A/V - IV) ffl 7ZV ffl TV (5) 

and the conjugation is generated by TV . 

In other words, the controlled map is 

U'\ ( 

- ' = e W e < v >etf> I , (6) 




where e {/} = e -{v} e Wv-rv} e {nv} e {rv}^ 
Proof: The conjugation equation 

e {ry} (Hmvmf) = Hmnv, 

has the solution 

/ = (-V) ffl (-77) me- {TV} (Hm KV). (7) 
In Appendix A we prove the identity 

e^W = U ffl Wm (-17), (8) 
for functions U,W E A. Applying Eq. (8) twice in Eq. (7) yields 

f = (-V) ffl (-if) ffl (-IV) ffl 77 ffl 7£V ffl TV 

= (-v) ffl (-e -{H} rv) ffl nv ffl rv, 

which reduces to Eq. (5) as e~^T = T — J\f by the definition of A/". □ 
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Note that the map generated by H EH TZV is not integrable in general (since 
TZV does depend on the angle tp). However, we point out two specific cases 
where it is integrable: The first one is when the frequency w/(2ir) is irrational 
and the second one is in the resonant case when V does not contain any mode 
k 7^ such that ktu e Z. In these two cases, TZV depends only on the action 
A. Therefore, the controlled map generated by H EH V EH / is integrable since 
it is canonically conjugate to H EH TZV = H + TZV . 

The Baker-Campbell-Hausdorff formula yields the expansion of the control 
term / : Its term linear in V which is —V + (A/V — TV) + TZV + TV, vanishes, 
so that its dominant term is of order V 2 : 

/ = -im^V + \{V}TZV + 0(V% (9) 

Note that in the continuous case, the dominant term of the control term 
was [13,15] 

TZ + 1 

By introducing the time step r between two iterations of the map, we define 
the operator 

1 _ e - T W 
Ti T = , 

r 

and T r as the pseudo-inverse of Ti T by Eq. (4). The projectors TZ T and Af T are 
defined accordingly The control term is given by 

f T = r- 1 [(-tV) EH (tKV - T T V) EH tTZ t V EH T T V] . 

The dominant terms of f T are given by 

fr = -{r T v}^^v + \{v}tz t v + o(v 3 ). 

Therefore when r tends to zero, the dominant term of the control term for 
symplectic maps tends to the one of the time-continuous Hamiltonian systems. 

Remark : If we choose the controlled map to be generated by H EH / EH V, 
then the control term / is defined from the control term / as / = e^f. 
The dominant term of / coincides with the one of / (higher order terms do 
not coincide in general). If the controlled map is chosen to be generated by 
/ EH H EH V, the control term / is / = e {H} f. 
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3 Maps in mixed coordinates 



We consider area-preserving maps obtained from a generating function of the 
form 

S(A', tp) = A'p + H(A') + eV(A', ip). 
In this section, we kept e for bookkeeping purposes. The map reads 



A = A' + ed v V(A',<p), (10) 
p' = ip + u(A')+ed A V(A',ip), (11) 

where oo(A) = d A H. Here 8aV(A', p) denotes the partial derivative of V with 
respect to the action (first variable) and d lfi V(A', p) denotes the partial deriva- 
tive of V with respect to the angle (second variable). 

Note that for e — 0, the map (10)-(11) is integrable and identical to the map 
generated by the Lie transform e^- H \ Thus we define the linear operators Ti., 
T, 1Z and Af as above. 

Our aim is to modify the generating function with a control term of order e 2 
such that the controlled map is closer to integrability than the original map. 
We consider the controlled generating function 

S C (A', ud) = A'p + H(A') + eV(A', ip) + e 2 f(A', ip). 

To construct the control term /, we perform a near-identity canonical trans- 
formation with generating function X(A , ip) = A ip + ex{Ao, ip) which maps 
(A,tp) into (A ,p ) and (A',tp') into (A' ,p' ) using the same function x- The 
scheme for the change of coordinates of the original map is depicted on Fig. 1. 
The canonical change of coordinates maps the system obtained by the gener- 
ating function S c into a system obtained by the generating function So : 



S (A' , ip ) = A' ip + H(A' ) + ez(A' , p ) 

d v z(A' , Vo )d A x{Ao, po) + d vX (A, + Lo(A' ))d A V(A' , Vo ) 

+ 0(e 3 ), (12) 



+e 2 



+\cu'(A' ) (d vX (4» + ui(A> ))) 2 + f(A' , ip ) 
where 

z(A' , ip ) = V(A' , ipo) + xK, ipo + wK)) - xK, <p ). (13) 
We refer to Appendix B for the detailed computations. 
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x 



X 



(4»<Po) 



(4,',<Po') 



Fig. 1. Diagram of the generating functions for the canonical changes of variables. 



Now we choose the generating function x such that the order e in Eq. (12) 
vanishes. To determine Xi we expand V as V = 1ZV + AfV. We require that : 

X (A, <p) - X (A, y? + cu(A)) = AfV (A, tp). 

By expanding V in Fourier series, i.e. V(A, <p) = Y,k&Vk(A)e lk,fi , this reads 

xM.*>)= E iS e * (14) 

fc S.t. ui{A)k0irZ 

which means that x(A, tp+u(A)) = -TV (A, <p). As x(A, ip+uj(A)) = e {H} x(A, ip), 
we find 

X = AfV - TV. (15) 

The control term is chosen such that the order e 2 of S given by Eq. (12) 
vanishes. Its expression 

f(A, <p) = (d A V) d v TV - ^ {d^TVf + (d v KV) d A {NV - TV), (16) 

represents the dominant term (of order e 2 ) of a control term that would cancel 
all the terms of higher orders in S . Using this control term, the controlled map 
obtained from the generating function S c is conjugate to the map generated 
by 

S (A', <p) = A'ip + H(A') + eRV + 0(e 3 ). 

As in the preceding section, in the two cases where u is irrational or V does not 
contain resonant modes, the map generated by So is e 3 -close to integrability. 
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1 2 3 4 5 6 

* 



Fig. 2. Phase portrait of the standard map S for e = 1.2. 
4 Numerical examples 



^.1 Application to the standard map 



The standard map S is 

A' = A + e sin </?, 
ip' — (p + A' mod 2?t. 

A phase portrait of this map for e = 1.2 is given in Fig. 2. There are no 
Kolmogorov-Arnold-Moser (KAM) tori (acting as barriers in phase space) at 
this value of e (and higher). The critical value of the parameter e for which 
all KAM tori are broken is e st d ~ 0.9716. 

The standard map is obtained from the generating function in mixed coordi- 
nates 

A' 2 

S(A', (p)=A'(p + — + e cos ip, 

i.e. V(A,(p) = ecosp and H(A) = A 2 /2. The generating function x given by 
Eq. (14) is thus 

, a \ sin (v? - A/2) 

The dominant control term given by Eq. (16) is 

2 cos% + A/2) 
f^,MM = -e 8sin2(A/2) , (17) 

and the resulting map iS m i Xj 2 generates the phase portrait displayed on Fig. 3. 
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Fig. 3. Phase portrait of the controlled standard map <S m i x ,2 with the control 
term (17) for e = 1.2. 

Note that besides the set of invariant tori which have been created in between 
the two primary resonances, the modification of phase space is significant. This 
comes from the fact that the control term / m i X; 2 does not induce a small mod- 
ification of the standard map : actually, it is unbounded when A approaches 
a primary resonance, i.e. when A G 2nZ. 

To obtain a smaller control term, we will localize it near a region where A is 
close to Ti modulo 2n. This corresponds to the centered resonance approxima- 
tion [24,25,19]. If one keeps only the leading term, the control term becomes 
independent of A (which may also be easier to implement in practice) : 



The controlled standard map <S c , a generated by f a is 

A' = A + e sin (p H sin 2(p : 

8 

<f' — ip + A' mod 2n. 

The phase portrait of <S Cia is displayed on Fig. 4. It clearly shows that adding 
the control term / a created a lot of invariant tori in the region near A = it 
whereas the structures near the primary islands are essentially preserved. 

If we consider higher order terms in the expansion near A = tt mod 2-7T, the 
control term becomes 



f b (A,<f) = sinV sin2v9sinA + cosVcos 2 (A/2)) + 0(cos 3 (A/2)). 



(18) 



1 




(19) 
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Fig. 4. Phase portrait of the controlled standard map 5 Ci(l with the approximate 
control term (18) for e = 1.2. 

If we neglect the cubic order, the controlled map iS Cj b becomes 



e 2 

A = A' - e sin ip (sin 2ip sin 2 ( A'/ 2) - cos 2(p sin A') , 

8 

e 2 

ip' = (p + A' + — (sin 2ip cos A 1 + cos 2 (p sin A'). 

The phase portrait of <S Ci b is very similar to the one of <S Cia shown on Fig. 4. In 
addition to the set of invariant tori created in the region near A = it mod 2ir 
in a very similar way as for <S C;a , it is worth noting that the invariant tori of S c ^ 
are more robust than the ones of <S C;a (i.e. they survive to slightly higher values 
of e). More precisely, iS Cja has invariant tori up to e a ~ 1.68, whereas invariant 
tori of iS Cj b survive up to ~ 1.81. Note that the controlled map <S m i Xi 2 has 
invariant tori up to £ m i x ,2 ~ 1-66 while the standard map has invariant tori 
up to £ std ~ 0.97. These values were obtained using Laskar's frequency map 
analysis [26,27,28,29]. 

In order to compare the control term obtained with mixed coordinates and 
the one using the Lie transform, note that the standard map is obtained by 
a Lie transform using H(A) = A 2 /2 and V(ip) = ecosy?. For A ^ 2ttZ, we 
notice that 1ZV = 0. The action of the operator r on V becomes 

sin(y9 + A/2) 
2sin(A/2) • 

The dominant order of the control term given by Eq. (9) becomes 

2 / \ 2 



Lie,2 



sin 



(A/2) J 
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Since the controlled map is difficult to explicit analytically, we consider the 
simplified case where we consider only the region near A = ir. The approxi- 
mation of the control term becomes 

e 2 

/Lie,a = -g" SmV- 

The controlled map is given by Eq. (6), which is 



A^A-dyVM-dyfubM, 

if' = if + lo(A - d v V((p) - d^fucM), 

since V and fu e ,a do not depend on the action A. The associated controlled 
map is thus the same as the one obtained with mixed coordinates <S Cja . 



4-2 Application to the tokamap 



The tokamap [21,22] has been proposed as a model map for toroidal chaotic 
magnetic fields. It describes the motion of field lines on the poloidal section 
in the toroidal geometry. This symplectic map {A, if) i— > [A', if'), where A is 
the toroidal flux and <f is the poloidal angle, is generated by the function 

S(A', f) = A'if + H(A') - e -^w{<p), 
with w(<f) = cos if. It reads 



A' 

A' = A — e — sin f, 

1 + A' r 

where q(A) = l/cu(A) is called the g-profile. In our computation, we choose 
u(A) = d A H = l/q(A) = tt(2 -A)(2-2A + A 2 )/2 and e = 9/(4tt). A phase 
portrait of this map is shown in Fig. 5. 

The controlled map must satisfy some specific requirements : In addition to 
being area-preserving, the toroidal geometry must be preserved which means 
that A has to be positive for all times (for a circular section, A = r 2 where r 
is the dimensionless radial coordinate). 

The generating function determined by Eq. (14) is 

eA sin(y? -u(A)/2) 



2(1 + A) sm(cu(A)/2) 
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Fig. 5. Phase portrait of the tokamap for e = 9/ (47r). 
The control term given by Eq. (16) is 



e 2 A cos(y + ^(A)/2) 
71 2(1 + A) 2 sin(w(A)/2) 



cos<^ cos(</? + u(A)/2) 

TTA ~ 4sin(a,(A)/2) 



We simplify the control term by considering the region near A = 1/2. To 
preserve positive values for A, we keep the term A/(l + A) in front of the 
control term. The control term becomes 

„ , , \ e 2 A cos(</? + a) (2 ll7rcos(<z> + a)\ . . 

/ a (A, tp) = -—J^T " ( o COS V 9 + • " , 20 

l + i4 3sma \3 64 sma J 

where a = 157r/32. It means that the potential w of the tokamap is now 
slightly modified into 

, . cos(y? + a) (2 117T cos(v? + a)\ 

w a ((p) = cosif + e—-. -cosv? + — : • 

3 sma \3 64 sma J 

The functions w and w a are plotted on Fig. 6. A phase portrait of the tokamap 
for e = 9/(47r) with the control term given by Eq. (20) is shown on Fig. 7. We 
clearly see that a lot of invariant tori have been created by the addition of the 
control term in the region A = 1/2. This shows that a small and appropriate 
modification of the map (through the potential function) can drastically reduce 
the chaotic zones and hence the transport properties of the system. 
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A Proof of e^w = U ffl W ffl (-U) 



First we recall the identity e^j3 = e a (3e~ a for linear operators a and /3, where 
{•} is the commutator, i.e. {a} (3 = af3 — f3a. We apply it to a — {U} and 
(5 = {W} in L(A), where U and W are functions in A. Then {{f/}} is the 
linear operator acting on L(A) as the commutation with the element {U}, i.e. 
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{{U}}{W} = {U}{W} - {W}{U}. Thus 

e {{u}} {w} = e {u} {w}e -{u} 

On the other hand, we have 

eMe {U} {W}e~W) = e^V^e"^ = e^^ w ^~ u \ 

where the second equality follows from definition of the warped addition © 
given by Eq. (2). 

The Jacobi identity implies {{U}W} = {U}{W} - {W}{U} = {{U}}{W}, 
from which it follows by recursion that {{£/} T W} = {{^ }} n {W}, for all 
n G Z, and hence 

e« u »{W} = {e%}. 
Since we have {U} © {W} = {U ffl W} by definition (3), we obtain Eq. (8). 



B Expression of the generating function S 



In what follows, the functions Va and denote the partial derivatives of any 
function V with respect to A and (p. In the same way, Va^ denotes the second 
partial derivative OaO^V. The determination of So follows the scheme depicted 
on Fig. 1. The computation is done mainly in two steps : First, we derive the 
generating function which maps the variables (A, ip) to (A' , ip' Q ) by eliminating 
the dependency on the variables (A', <//). Secondly, we derive the expression 
of the generating function Sq which maps the variables (A , (p ) to (A' , (p' ) by 
eliminating the dependency on the variables (A, if). 

The variables (A' , ip' ) are defined implicitly as functions of (A', ip') by the 
canonical transformation generated by x '■ 



A' = A' + e Xv> (Aoi<P')i (B.l) 
<p' Q = <p' + e X A(A! Q ,ip'). (B.2) 

The generating function S c gives 



A = A' + eV v {A', <p) + e 2 f v (A', v ), (B.3) 
^ = p + uj(A') + eV A (A', ip) + e 2 f A (A', ip). (B.4) 

First, we rewrite Eq. (B.4) by using Eq. (B.l) in order to get the expression 
of if' in terms of the variables (A' , if) up to order e 2 : 
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tp' = if + u(A ) + e [u A {A! Q ) X ^ P + wK)) + VkK, </?)] 



+5 Z 



1 



We substitute the expression of <// into Eq. (B.2) : 



0(e 3 ). 



<p' = tp + u(A' ) + e [u A (A' ) Xip (Ao> V? + "(4>)) + XaK, <p + w(^)) + Va(^o, </?)] 



uja{A' ) < 2 Xv {A' , p + wK))x w K, y? + 



+V AA (A' , <p)x v (4>, V + w(^)) + VkK, v?)xa^(K, + 
+w A (4>)x v K, V? + p + ^K)) + /aK, + 0(e 3 ). 



The substitution of the expression of tp' into Eq. (B.l) gives A' as a function 
of and tp. Inserting this expression into Eq. (B.3) gives 



A = A' + e [ X<P (A' , p + lu(A' )) + V^(A' , tp)} 

+e 2 [V Aip {A' , p)x<p(A'o, V? + w(^)) + VkK, y?)x w K, + ^K)) 

y? + ^K)) Xw K, y? + uj(A' )) + /^K, v?)] + 0(e 3 ). 

The expressions of A and ip' as functions of A' and </? can be obtained by the 
generating function : 



S (A' , p) = A' tp + H(A> ) + s [V(A' , tp) + X (A' , p + u(A' ))\ 



V A (A' ,p)x v (A' ,tp + u(A' )) 

\0(e 3 ). 

Next, we eliminate the variables A and p using the equations : 



+ -lu a (A' ) X 1(A' , p + u;(A' )) + /K, y?) 



A = A + ex v (A ,p), 
p = p + ex A (A ,p). 



(B.5) 
(B.6) 



In order to perform the same type of procedure to eliminate A and tp, we first 
need to invert Eqs. (B.5)-(B.6) : 
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A = A - ex v (A, ip Q ) 

+e 2 \x v (A, ¥o)xa^A, (p ) + Xa(A, </?o)x w (A ip )} + 0(e 3 ), 
tp = tp - exa(A, (p ) 

+e 2 [ XA (A, ¥o)xa v (A, <Po) + X<e( A , <Po)xaa( a , ip )} + 0(e 3 ), 

which can also be obtained from the generating function 

X(A, v? ) = A<p - e X {A, Vo ) + e 2 XA {A, (p )x v (A, Vo) + 0(e 3 ). 

The same type of substitutions as for eliminating (A', ip') leads to the expres- 
sion of A and ip' as functions of A' and (po : 

A = A' + e \x<p(A' , y? + ^K)) - X<p( A 'o, <Po) + K( A 'o, <Po)] 

+e 2 V Aip (A' , <p Q )x v (K, V?o + uj{A' )) + V A (A' , Vo ) Xw {A' , y? + u>(A' )) 

-V W (A' , ip )x A (A' , <p ) - ip + uj(A' ))xa(A' , <p ) 

+u A (A' )x< p ( A 'o, V?o + iv(A' )) X<pip (A' , V?o + w(^)) 
+Xa v (^o, ¥>o)(x v (4,, ¥>o) - X^o, V?o + wK))) 

(f/ = (p + uj(A' ) + e [V A (A' , v^o) + c^K) X ^K, V?o + ^K)) 
V?o + ^K)) - XaK, <A))] 

+e 2 Xa(^o, ¥>o)xa v (^,, Vo) + X<p(A' , tPo)XAA( A ' , <Po) 

-Xa( a ' , (Po)V Aip (A' , v? ) - uj a {A' ) Xa {A' , (p )xw( A 'o, ¥o + u(A' )) 
-Xa( a ' , <Po)XA<p( A ' , + ^K)) + X^K, + uj(A' ))V aa (A' , y? ) 
+V A (A' , <p ) (XA<p(A'oi V?o + wK)) + w^^XwK, V?o + w(^))) 
+u; A K) X ^K, v? + ^K)) X a^K, V?o + w(^)) 
+u; A K) 2 X^K, + ^K)) Xw K, y? + 



+-u AA (A' )x 2 JA / , Vo + u(A' )) + f A (A' , v9 ) 



+ 0(e 3 ). 



The expressions of A and ip' can be obtained using the generating function 



S (A' , v? ) = A' oVo + if (4,) + e [V(A' , Vo ) + X (A' , y? + uj(A' )) - X ( A ' , <Po)] 



(V A (A' , ip ) - X a(A' , <p )) X<p ( A o, V?o + uj(A' )) 



~ (V<p( A 'o, <Po) ~ xA A i <A))) Xa(^o, y?o) 
+lu A (A' ) x 2 JA' , v9 + c^K)) + f(A' Q , <p ) 



+ 0(e 3 ). 



Using Eq. (13) yields Eq. (12). 
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